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We use probabilistic methods to study classical solutions for sys- 
tems of interacting semilinear parabolic partial differential equations. 
In a modeling framework for a financial market with interacting Ito 
and point processes, such PDEs are shown to provide a natural de- 
scription for the solution of hedging and valuation problems for con- 
tingent claims with a recursive payoff structure. 

1. Introduction. Reaction-diffusion systems are systems of semilinear 
parabolic partial differential equations which can interact in a possibly non- 
linear way. They appear as models for phenomena from various areas of 
applications, ranging from ecological systems and biological pattern forma- 
tion to chemical reactions; see Smoller (1994) for references. This article is 
concerned with applications to hedging and valuation problems in math- 
ematical finance. Standard existence and uniqueness results for reaction- 
diffusion systems may not apply here because the coefficient functions of 
typical parametrizations in finance often are unbounded or do not satisfy 
linear growth constraints. A first contribution of this article is to address 
that issue by proving results on classical solutions in a fairly general context. 
A second contribution is an integrated treatment of contingent claims in the 
context of a flexible Markovian framework which incorporates new features 
and includes more specific models studied so far. We use the developed PDE 
techniques to provide results on the valuation and hedging of claims with a 
recursive payoff structure. 
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The article is structured as follows. Section 2 contributes existence and 
uniqueness results for classical solutions of reaction-diffusion systems by 
showing existence of a unique fixed point of a suitable Feynman-Kac rep- 
resentation, provided that the underlying diffusion stays within a given do- 
main and the coefficient functions satisfy local conditions. This first yields 
results under a (global) Lipschitz condition on the interaction term, and 
these are extended to a local Lipschitz condition by exploiting an additional 
monotonicity assumption on the interaction. The latter arises naturally in 
valuation problems from mathematical finance. Section 3 introduces for the 
subsequent applications a stochastic model with interacting Ito and point 
processes and gives a construction by a change of measure. The model con- 
sists of a system of stochastic differential equations which describes the 
Markovian dynamics of an Ito process S and a further finite-state pro- 
cess rj driven by the point processes. This SDE system is nonstandard in 
that the driving process can itself depend on the solution, similarly as in 
Jacod and Protter (1982). 

In Section 4, this framework is used as a model for an incomplete financial 
market, with S describing the prices of tradable assets, for example, stock 
indices, while the process rj represents further (not directly tradable) sources 
of financial risk, for example, rating and credit events or the state of an in- 
surance contract. An important feature is that our model allows for a mutual 
dependence between S and r], in that the drift and volatility of S can depend 
on the finite-state process rj while the intensities for changes of rj can in turn 
depend on the current value of S. In the context of mathematical finance, 
this can be seen as both a fusion and a generalization of a Markov chain 
modulated diffusion model of Black-Scholes type, as in Di Masi, Kabanov 
and Runggaldier (1994), and of the Cox process model for credit risk from 
Lando (1998) or the conditional Markov chain model, respectively. Another 
contribution is that we not only study a pure pricing approach under an a 
priori given pricing measure, but use a combination of valuation and hedging 
ideas to determine both a valuation and a locally risk-minimizing hedging 
strategy. Following Duffie, Schroder and Skiadas (1996), we allow for claims 
whose payoffs can depend not just on the state variables S and t] but also on 
the valuation process of the claim itself. This leads to a recursive valuation 
problem, and it turns out that reaction-diffusion systems as in Section 2 
provide a natural, convenient and constructive description in terms of PDEs 
for the solution to valuing and hedging problems for such claims. Section 5 
discusses application examples and possible extensions, mainly with a view 
towards hedging of credit risk. This illustrates the flexibility of our model 
and results and also relates them to the existing literature on this topic 
which gained much interest recently; see Jeanblanc and Rutkowski (2003) 
for an overview and more references. 
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2. Classical solutions for reaction diffusion systems. In this section, we 
use stochastic methods to derive existence and uniqueness results for classi- 
cal solutions of interacting systems of semilinear parabolic partial differential 
equations (PDEs). Such systems are also known as reaction-diffusion equa- 
tions and play in subsequent sections a key role in our solutions to various 
valuation and hedging problems from mathematical finance. There we con- 
sider a Markovian setting where an Ito process S models the prices of the 
tradable assets, and further nontraded factors of risk are represented by a 
finite-state process r/. Similarly as in the Black-Scholes model, the solutions 
to our valuation and hedging problems can be conveniently described via 
PDEs. But the nontradable factors lead to an interacting system of PDEs; 
each single PDE corresponds to a possible state of r], and the interaction 
between the PDEs reflects the impact from the evolution of i] on the valua- 
tion. 

We first derive results for PDE systems where the interaction satisfies a 
global Lipschitz condition, and extend these to a type of monotonic local 
Lipschitz interaction. The latter is relevant for the applications to valuation 
and hedging problems. We strive for general assumptions on the coefficient 
functions which are satisfied by typical financial models. 

2.1. General framework. Fix m E N, a time horizon T € (0,oo) and a 
domain (open connected subset) D in M.'^. For each starting point {t,x,k) G 
[0, T] X D X {1, . . . , m}, consider the following stochastic differential equation 
(SDE) in R'^: 

(2.1) 

^xt,.,k ^ r,(,, x*'-''=) ds + j2 ^kA^^ ^^'"'') dwi, s G [t, T], 

i=i 

for continuous functions Ffc :[0,T]xD-^ and T^kj :[0,T] x D -^R'^, j = 
1, . . . ,r, with an M^-valued Brownian motion W = (W^)j=i^,,,^r- We write F^ 
and each j- as a dxl column vector and define the matrix-valued function 
Efc :[0,T]xD^ M'^x^- by ^ := (Sfcj)\ For any k, Ffe and S^,,, j = 1, . . . , r, 
are assumed locally Lipschitz-continuous in x, uniformly in t: 

For each compact subset JC oi D, 

there is a constant c = c(/C) < cx3 such that 

(2.2) 

\G{t,x)-G{t,y)\<c\x-y\ 

for all t G [0,T],x,y e K and G G {Ffc,Sfc,i, . . . 

By Theorem V.38 in Protter (2004), condition (2.2) implies that (2.1) has a 
unique (strong) solution for any given tuple (fi, .7^, F, P, VF) up to a possibly 
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finite random explosion time. We impose the additional global and proba- 
bilistic assumption that for all {t,x,k), the solution X*'^''' does not leave D 
before T, that is, 

(2.3) P[X*'^'*^ € D for all s G [t, T]]=l. 

This includes that X*'^'*^ does not explode to infinity: 



sup IX*'""'*-'! < oo 

s&[t,T] 



1. 



By Theorem II. 5. 2 of Kunita (1984), (2.2) and (2.3) imply that X*'^''' has a 
version such that 

(2.4) {t,x,s) 1-^ X*'^'''' is P-a.s. continuous. 

2.2. Fixed points of the Feynman-Kac representation {generalized solu- 
tions). For existence and uniqueness problems of nonlinear PDEs, it is 
common to consider generalized solutions, namely solutions of a correspond- 
ing integral equation. These in general require and possess less regularity, 
and additional assumptions are needed to ensure that a solution to the in- 
tegral equation is also a classical solution to the PDE. See Chapter 6.1 of 
Pazy (1983) for an analytic version of this approach and Preidlin (1985) for 
a probabilistic version. 

For the PDE (2.8) that we consider below, the integral form of the stochas- 
tic approach is the well-known Feynman-Kac representation. Since the PDE 
is nonlinear, the solution itself appears within the expectation so that we 
have to look for a fixed point. To make this precise, we start with continu- 
ous functions h: D ^W", g ■.[0,T] x D xW^ and c:[0,T] x D ^ R"'. 
Given sufficient integrability, one can then define an operator F on functions 
V by 



(Fv)''it,x) :=E 

(2.5) 



+ 



T 
T 



with k = 1, . . . ,m, (t,x) G [0,T] x D. Under suitable conditions, F has a 
unique fixed point: 



Proposition 2.1. Assume (2.2) and (2.3) hold. Let h, g and c be con- 
tinuous, with h and g bounded, and with c bounded from above in all coor- 
dinates. Then F defines a mapping from the Banach space Cf, := Cb([0,T] x 



REACTION-DIFFUSION SYSTEMS AND HEDGING 5 

D,M"^) of continuous bounded functions v:[0,T] x D ^ into itself. As- 
sume further that {t, x, v) i— > g{t, x, v) is Lip schitz- continuous in v, uniformly 
in t and x, that is, there exists L < oo such that 

\g{t,x,vi) - g{t,x,V2)\ < L\vi - V2\ 

(2.6) 

for all t € [0,T], x e D and vi,V2eW^. 
Then F is a contraction on Ch with respect to the norm 
(2.7) \\v\\p-= sup e-^^^-^\v{t,x)\ 

(i,x)G[0,T]xD 

for 13 < oo large enough. In particular, F has then a unique fixed point v € 

Proof. By the boundedness assumptions and (2.3), (2.4), the operator 
F maps Cb into itself. The norm (2.7) is equivalent to the usual sup-norm 
on Cb, and by assumption, c has values in {—oo,K]"^ for some constant 
K G [0, oo). Using this and (2.6), we obtain for v,w &Cb and (3>0 that 

e-/3{7^-t)|(Fz;)^(t,x) - {Fw)''{t,x)\ 
E 



1 



cP(T-t) 



pKT 

< F 

- e/3(T-t) ^ 



(/(s,X*•-'^^;(s,X*'=^'*^)) 

- /(s, «;(s, X*'"''=)))e/t' -'(«'^"^'')'^- ds 
\g\s,Xl^-^\v{s,Xl^-^^)) 



g''{s, X'f''',w{s, x*'^.fc)) |e-/3(r-'*)e/3(^-^) ds 

KT rT 



1/3 



< 7, — \\V — W\ 

(3 



/3 



for all (t,x) G [0,T] x D and k = 1, . . . ,m. Thus i*" is a contraction for (3 > 

□ 

2.3. Classical solutions under local regularity. Now define operators C^, 
k = 1, . . . ,m, on sufficiently smooth functions / : [0, T] x Z) — M by 

iC'f){t,x)=Y^nit,x)^it,x) + - Y: ai{t,x)^-l-{t,x) 

i=l *J=1 
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with 

Consider the following system of semilinear PDEs with k = l,...,m and 
boundary conditions at terminal time T: 

d 

— f '^(t, x) + C^v^{t, x) + c^{t, x)v^{t, x) + g^{t, X, v{t, x)) = 

(2.8) for € [0,r) X D, 

v''{T, x) = h^{x) for X G L>. 

These m PDEs are interacting via the g-teixn which may depend on all 
components of v{t,x) = {v^{t,x))k=i,...,m- Our goal is to show that the fixed 
point V from Proposition 2.1 is the unique bounded classical solution to (2.8). 
To this end we apply a Feynman-Kac type result from Heath and Schweizer 
(2001) that relies on classical results by Friedman (1975) and requires only 
local assumptions on the coefficient functions. Note that the subsequent 
results do not require further boundary conditions for the PDE (2.8); this 
is due to condition (2.3). We impose the following additional conditions on 
the coefficients of the SDE (2.1) and the PDE (2.8): 

There exists a sequence {Dn)n£N of bounded domains with 

(2.9) closure Dn ^ D such that UJJLi Dn = D, each D„ has a C^- 
boundary, and for each n and k= 1, . . . , m, 



(2.10) 



the functions and = S^.S^'' are uniformly Lipschitz- 
continuous on [0, T] x Z)„, 



(2.11) detafc(t,x) /O for all {t,x) G [0,r] x D, 

(2.12) {t,x) 1-^ c{t,x) is uniformly Holder-continuous on [0,T] x Z)„, 



(2.13) 



{t,x,v) I— > g{t,x,v) is uniformly Holder-continuous 
on [0,r] X Z)„ X M™. 



Remark 2.2. We aim for a classical solution, not a generalized solution 
(as already provided by Proposition 2.1) or a solution in some other weak 
sense; see, for example, Pardoux (1999) for viscosity solutions of similar PDE 
systems. To the best of our knowledge, the subsequent results on classical 
solutions have not been available so far under the assumptions given here. 

Apart from the global condition (2.3), we only impose mild local con- 
ditions on the coefficient functions; we assume them neither bounded nor 
satisfying a global (linear) growth condition over the possibly unbounded do- 
main D. This is crucial because such restrictive assumptions could exclude 
parametrizations which are typical in financial models; the only global con- 
dition (2.3) is probabilistic and must be verified on a case by case basis. 
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for instance by means of Feller's test for explosion. An example is given in 
Heath and Schweizer (2001). 

Let us denote by := Cjj ([0,r) x L>,M™) the spaces of continuous 
(bounded) functions : [0, T] x D —>■ R™ which are of class C^'"^ with respect 
to (t, x) £ [0, T) X D. Note that the C-'^'^-condition is imposed only on [0, T) x 
D while continuity is required on all of [0, T] x D. 

Proposition 2.3. Assume that (2.9)-(2.13) hold in addition to all the 
assumptions for Proposition 2.1. Then the system (2.8) of interacting semi- 

1 2 

linear PDEs has a unique classical solution in C^' , which is given by the 
fixed point v from Proposition 2.1. 

Proof. Recalling that is a contraction on Cb by Proposition 2.1, 
we first show that F maps bounded functions w that are locally Holder- 

1 2 

continuous in {t,x) on [0,T) x D into C^' , and that for each such w, the 
function v := Fw satisfies the following system of m PDEs with terminal 
conditions: 

d 

-K-v^{t, x) + d'v^{t, x) + c^{t, x)v''{t, x) + g^{t, X, w{t, x)) = 0, 

(2.14) 

{t,x) G [0,r) X D, 

v^{T,x) = h^[x), x(£D. 

It is evident from the definitions of v and F that v satisfies the terminal 
condition and is bounded, with a bound on v that depends only on the 
bounds for h, g and c. To prove the above two assertions, it suffices to 
show for any e > that v is in C^''^{[0,T - e) x D,W^) and satisfies (2.14) 
on [0,T -e) X D instead of [0,r) x D. So fix arbitrary e G (0,r) and k G 
{1, . . . , m}, and let T' ■.= T — e. For any (t, x) G [0, T'] x D, conditioning on 
Tt' gives 

v''{t,x) = {Fwf{t,x) 



E 



E 



+ 



E 



+ 



Jt 
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by using the Markov property of X^'^'^ for the last equahty; compare the 
proof of Theorem 1 in Heath and Schweizer (2001). Using that Theorem 1 
and the above representation for v'^{t,x) on [0,r'] x D, we obtain that 
{v^)k=i,...,m is in C^'2([0,r') x D,W^) and satisfies the PDE 

^v'^it, x) + C^v'^it, x) + c^it, x)v^{t, x) + g^it, x, w{t, x)) = 0, 

(2.15) * 

{t,x) € [0,T') X L>, 

if we can verify the assumptions [A1]-[A3] from Heath and Schweizer (2001). 
[Al] and [A2] are precisely (2.2) and (2.3). We check the list [A3'] of condi- 
tions in Heath and Schweizer (2001) whose combination implies [A3]. Condi- 
tions [A3'], [A3a'] and [A3c'] in Heath and Schweizer (2001) are exactly (2.9), 
(2.10) and (2.12). By Lemma 3 in Heath and Schweizer (2001), the conti- 
nuity of Sfc in combination with (2.11) implies their condition [A3b'] (a is 
uniform elliptic on [0,T] x !)„), and v £ implies [A3e'] {v is finite and 
continuous). To verify [A3d'] [i.e., (2.16) below], note that w is by assump- 
tion uniformly Holder-continuous on the compact subsets [0,r'] x Hence 
(2.13) implies that the composition 

(2.16) (t,x) g{t,x,w{t,x)) is uniformly Holder-continuous on [0,r'] x l)„ 

which is [A3d']. Since e > was arbitrary, we conclude from (2.15) that 
V = {v'')k=i,...,m = Fw is in Cl'^ and satisfies the PDE (2.14) on [0,r) x D. 

It follows that the fixed point w € Cf, from Proposition 2.1 is approximated 
in Cf, (i.e., in the sup-norm) by a sequence (fn)nGNo = {Fvn-i)neN from C^'^ 
if we choose vo = w locally Holder-continuous in {t,x). To prove that v is an 
element of C^' C Ch and satisfies (2.8), it suffices by the preceding argument 
to show that 

(2.17) -0 is locally Holder-continuous in {t,x) on [0,T) x D. 

To establish (2.17), we employ an a priori Holder estimate for the sequence 
{vn)n£'N which is local in (t,x) but uniform in n. Let Q denote a bounded 
domain with Q ^ [0,r) x D, and let Q' be some sub domain of Q having a 
strictly positive distance to dQ fl (0,T), where dQ denotes the boundary of 
Q. Then there is a Q'-local Holder estimate which holds uniformly for all 
functions of the sequence (f^). More precisely, there exist by Theorems 6 
and 7 in Chapter 4.2 of Krylov (1987) some constants f3 € (0, 1) and N < oo 
such that the estimate 

\v':{t,x)-v'^{t',x')\ 

(2.18) 

< A^(lb^llL«>(Q) + ||/(-, t^n-l(-)) IIl^+1(Q))(|x - X'\ + \t - t'\'/'f 



REACTION-DIFFUSION SYSTEMS AND HEDGING 



9 



holds for all {t,x), {t',x') € Q', n G N and k = l,...,m. This uses that (2.14) 
holds for V = Vn = Fvn-i and w = Vn-i', note also that C^'^((5) is contained 
in the Sobolev space 

W^i+i(Q) of Krylov (1987) because we can choose Q 
to have smooth boundary and use approximations by C°°-functions. By the 
boundedness of Q and the and L'^+^-norms of the functions v^{t,x) 
and g''{t,x,Vn-i{t,x)), respectively, with respect to Lebesgue measure are 
bounded uniformly in n G N. In fact, we have (as already noted) a uniform 
bound on all |j^^.J^||L°°(Q) terms of h, g and c, and 

||/(-,^^n-i(-))llL<i+i(Q) < I|1|Il''+i(Q)II/IIl-(Q) <oo for k = l,...,m. 

Hence the first bracket on the right-hand side of (2.18) is bounded uniformly 
in n. For any such Q' , this yields a Q'-local Holder estimate for the sequence 
{vn)ne'N that is uniform in n. It follows that the (uniform) limit v Cb is 
locally Holder-continuous. This establishes (2.17). 

Uniqueness follows by the usual Feynman-Kac argument. In fact, we 
apply Ito's formula to the process v'' [s, Xj.'^''^) ex.p{J^ {u, X^^''') du) , s G 

[t,T), and use the PDE to show that any solution v G C^'^ to (2.8) is given 
by the Feynman-Kac representation and therefore a fixed point of F. Since 
the fixed point is unique, this yields v = v. □ 

2.4. Classical solutions under monotonic interaction. For our applica- 
tions in the context of hedging and pricing in finance, we need a solution to 
the PDE system (2.8) with a function g G C^{[0,T] x D x ]R™,M™) which 
is usually unbounded on [0,T] x D x M™, but locally bounded in v in the 
sense that 

g{t,x,v) is bounded on [0,T] x D x IC 

(2.19) 

for any compact subset /C of M™, 
and satisfies a monotonicity assumption of the following type: 

There exist Ki,K2 G [0, oo) such that for all = 1, . . . , m we have 

g''it,x,v) < +Ki + K2\v\ 

(2.20) for ah t e [0,T], x e D, v e {w e R"'\w'' >w\ Vj / k}, 

g^{t,x,v) > -Ki-K2\v\ 

for ah t G [0,r], xeD,ve{we M."^\w'' < w^, Vj / k}. 

Thus we have an upper (or lower) bound on coordinate k of g{t,x,v), at 
most linear in \v\, if the argument v has its largest (or smallest) coordinate 
for index k. We further suppose that 

{t,x,v) 1-^ g{t,x,v) is locally Lipschitz-continuous in v, 

(2.21) 

uniformly in [t,x). 
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Typical examples from mathematical finance which satisfy (2.20) and (2.21) 
are 

m 

gHt,x,v) = 6''it,x,v) + J2x''^it,x){v^ -v' + f^it,x,v)), 

i=i 

(2.22) 

k = I, . . . ,m, 

as in Section 4, or [as in Becherer (2004)] 

g''it,x,v) = 6^ it,x,v) + J2 > 2;) - (e"("' -^'+^'' - 1) , 

(2.23) 

fe = 1, . . . , m, 

with functions A^^^' G C^n[0,2"] x -D,[0,oo)) and 6^j''^ G (^{[([O.^] x D x 
M™,M) for k,j = l,...,m, which are locally Lipschitz-continuous in v, uni- 
formly in {t,x), and with a > 0. In Section 4, we shall work with (2.22) and 
replace boundedness of 6^ and f^^ by a linear growth condition in v; this is 
still covered by (2.20). 

Under the above assumptions on g, we cannot apply Proposition 2.3 di- 
rectly since g is not bounded, is not (globally) Lipschitz-continuous in v 
and does not satisfy (2.13) in general. But we can still obtain the following 
result: 



Theorem 2.4. Suppose h and c are continuous functions with h bounded 
and with c bounded from above in all coordinates. Assume that (2.2), (2.3) 
and (2.9)-(2.12) hold. Suppose g is in and satisfies the local boundedness 
condition (2.19), the monotonicity condition (2.20) and the local Lipschitz 
condition (2.21). Then the system (2.8) of PDEs has a unique classical so- 



lution V 
tation 



1,2 



([0,r)x A 



and V satisfies the Feynman-Kac represen- 



\t,x)=E 



(2.24) 



+ 



= (s,X*'^'^^)(s,X*'^'*^))e/t 



ds 



for k = 1, . . . ,m and {t, x) G [0, T] x D. 



Proof. It suffices to prove the assertion for the case where c is nonpos- 
itive in all coordinates because the general case can be reduced to this by 
passing to the transformed function {t,x,k) 1— > exp(—K(T — t))v(t,x,k) for 
a suitable constant K G [0, 00). 
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Since h is bounded, there is a constant e [0, cxd) such that \h^{x)\ < 
K3 for all X and k. For this constant and Ki^2 from (2.20) with K2 taken 
with respect to the max-norm on M™, we first define a truncation-boundary 
function k : [0, T] — > [0, 00) by 

(K3 + Ki{T-t), wheni^2 = 0, 

I ^2 
and then ^ : [0, T] x D x — > by truncating the third argument: 

g{t, X, v) := g{t, x, ( min ( max {v'', -K{t)), +K(i)))fc=i,...,m)- 

Then g is bounded by (2.19), and Lipschitz-continuous on [0,T] x Z)„ x R"^ 
for every n since 5 is in C^. By (2.21) the function g is also Lipschitz in 
V, uniformly for {t,x) E [0,T] x D, and so Proposition 2.3 yields a unique 
bounded solution v for the PDE (2.8) with g instead of g. Moreover, v is the 
fixed point of F defined with g instead of g in (2.5). We show below that 

(2.25) \v''{t,x)\<K{t) for (t,x) G [0,r] X L> and A; = l,...,m. 

Admitting this result for the moment, we get g{t,x,v{t,x)) = g{t,x,v(t,x)) 
for all {t,x) G [0,r] x D by the definition of g. Hence, v also solves the 
PDE (2.8) with g instead of g and satisfies (2.24). To see that v is the 
unique bounded solution to (2.8), let w denote another bounded solution. 
By taking larger if necessary, we can assume that \w^{t,x) \ < K3 < K{t) 
for all k,t,x. Then both v and w solve (2.8) not only with g but also with 
g, and this implies v = w hy the uniqueness in Proposition 2.3 applied for g. 

To finish the proof, it remains to establish (2.25). Fix arbitrary {t,x,k) G 
[0, T] X D x {1, . . . , m} and define the stopping time 

r:=inf{sG [t,r]|{;'=(s,X*'^'^) <k(s)} AT. 

Then 

{)'=(s,X*'^'^(u;)) > k{s) for all (w, s) G p,r[[ 

while v^{t,X!i.'^'^) < k{t); in fact, we have equality for t <T and inequality 
for T = T since v^(T, •) = h^{-) < K^. Hence, the definition of g and property 
(2.20) of g imply 

(2.26) ~g''{s,Xl'-'''{Lo),v{s,Xl'^'''{u;)))<K^+K2K{s) for (.;,s)G[t,r[ 
and therefore by using c < and integrating k 



v^{t,x) = E 



E 



12 



D. BECHERER AND M. SCHWEIZER 



E 



ds 



<E 



k{t)+ j\Ki + K2^{s))ds 



= E[K{T) + {K{t)-K{T))] 
= K{t), 

where the second equahty uses the strong Markov property of X*'^'*^; see 
Heath and Schweizer (2001). This gives the upper bound in (2.25), and the 
lower bound is proved in the same way by using r := inf{s G [t, T]\v^{s, X\'^'^) > 
— k(s)} a T. This completes the proof. □ 



3. A model with interacting Ito and point processes. In this section, 
we introduce a flexible Markovian model for an incomplete financial mar- 
ket, give a rigorous construction for it and provide some further properties. 
This is used in the next section to derive explicit and constructive results 
for various hedging and valuation approaches in terms of reaction-diffusion 
systems. The entire model is given by a system of stochastic differential 
equations (SDEs) of two types. The price process S of the tradable risky 
assets is modeled by an Ito process. In addition, there are nontradable fac- 
tors of uncertainty and risk which are represented by a finite-state process 
r] driven by a point process. A distinctive feature of our model is that it 
permits mutual dependences between S and 77. More precisely, both S and 
r] enter the coefficients of the SDE for the dynamics of S", and at the same 
time, the intensities controlling the jumps of rj depend on the current value 
of S. 



3.1. Model setup and assumptions. All modeling in the sequel takes place 
on some filtered probability space (fi, F, P) with a filtration F = {J^t)te[G,T] 
satisfying the usual conditions and a trivial u-field J-q. All semimartingales 
are taken to have right-continuous paths with left limits. 

We start with m G N and a domain D in M'^ satisfying (2.9), for example, 
D = or D = (0,00)'^. Let [3,7]) be a solution of the following system of 
SDEs with values in D x {1, . . . 

(3.1) 5o€l), dSt = T{t,St,m-)dt + ^{t,St,iit-)dWt, 

m. 

(3.2) 7?o € {1, ... , m}, dr^t = E 0' " k)I{k}{m-)dN^' , 

k,j=i 
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where T : [0, T] x L> x {1, . . . , m} ^ M*^ and S : [0, T] x D x {1, . . . , m} ^ M'^'^'' 
are with respect to {t,x) G [0,r] x D, W = {W'')i=i^,„^r is an M''-valued 
(P, F)-Brownian motion and N = (iV'^-')fcj=i ,,,^„^ is a multivariate F-adapted 
point process such that 

(3.3) (iVf^) has (P,F)-intensity X''^{t, St) for k,j = 1, . . . ,m 

with bounded functions X''^ '■[OjT] x D ^ [0, oo). Note that the process 
counting the jumps of rj from state k to j is not N^^ , but / I|fc|(r;_) dN^K 
If L> C (0,00)'^, one can rewrite (3.1) as a generahzed Black-Scholes model; 
denoting dS / S = {dS^ / S^)i=i^,,,^ci-, lit^Xjk) = diag(l/x*)j=i_,,,^dr(i, x, fc) and 
a{t, X, k) = diag(l/x*)j=i^,,,^rfi;(i, x, k), we have 

(3.4) —1 = j{t, St, vt-) dt + a{t, St, Vt-) dWf 

The model (3.1)-(3.3) is a nonstandard SDE system because of its depen- 
dence structure. The coefficients in the SDE (3.1) for S contain r], and the 
intensities in (3.3) of the point process N driving 77 depend in turn on S. 
We shall comment below on construction and properties. 

To apply the PDE results from Section 2, we need further regularity 
assumptions on the coefficients of the SDE (3.1). Because {t,x) 1— > G{t,x,k) 
is on [0,r] x D for G G {r,S} and any k, the mappings x 1-^ G{t,x,k) 
are locally Lipschitz-continuous in x, uniformly in t and k. This implies as 
in Section 2 that there is a unique strong solution X^'^'^ to the SDE 

(3.5) 

dxl'''''' = r{s,xl'''''',k)ds + ^{s,xl'''''',k)dWs, se[t,T], 

for any (t, x,k) G [0,T]x D x {1, . . . , m} up to a possibly finite random explo- 
sion time. As in Section 2.1 [for Ti^{t, x) := T{t, x, k) and Sfc(i, x) := S(t, x, k) 
there], we assume again that X^'^'^ does not leave D during [0,T]; that is, 
we suppose 

(3.6) G D for ah s G [t, T]] = 1 for any t, x, k. 

Remark 3.1. Intuitively, the SDEs (3.5) are related to (3.1) in the sense 
that 5 could be constructed successively from one jump time of 77 to the next 
by "pasting together" appropriate solutions to (3.5), using techniques similar 
to those known from the construction of finite-state Markov processes or Cox 
processes. But such a construction of {3,7]) becomes tedious in details, and 
the construction via a change of measure described in Section 3.3 appears 
more convenient. 
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3.2. Markov property and uniqueness in distribution. A standard way to 
show the Markov property is to prove uniqueness of a corresponding (time- 
inhomogeneous) martingale problem. We give here a direct argument which 
is similar in spirit. For a horizon T' £ [0,T] and a function h G Ch(D x 
{1, . . . , m}, M), we consider the PDE system 

d 

= vt(t,x,k) + T{t,x,k) gTad^v{t,x,k) + ^ ^ a^^ {t,x,k)v^i^j{t,x,k) 

(3.7) 

+ Y,^'''it,x){v{t,x,j)-v{t,x,k)), {t,x)e[0,T')xD, 
i=i 

for k € {1, . . . ,m}, with {a^-')ij=i^...^d = cl:= EE*'', and terminal conditions 
(3.8) v{T',x,k) = h{x,k), xeD. 

For brevity, we use subscripts for the partial derivatives of v. By Theo- 
rem 2.4 there is a unique bounded classical solution v G C^'^([0,r') x D x 
{1, . . . ,m},]R) to (3.7) and (3.8) for any given T' and h. The essential mar- 
tingale argument for Proposition 3.3 is 

Lemma 3.2. For v given as above, the process v{t, St,i]t), t S [0, T'], is 
a martingale. 

The proof is mainly an application of Ito's formula and given in the Ap- 
pendix. In the same way, another application of Ito's formula yields that 
{S, rf) solves the following martingale problem: For any continuous function 
/(x, k) on D X {1, ... , m} with compact support that is of class in x, the 
process 

f{Sum)- f{So.m)- f' Asf{Ss,vs)ds, te[o,r], 

Jo 

is a martingale, with the operators As being given by 

d 

Asf{x,k)=r{s,x,k)gTad^f{x,k) + ^ ^ a^^ {s,x,k)f^r^j{x,k) 

i,i=l 

(3.9) 

TO 

+ Y^X>'^{s,x){f{x,j)-f{x,k)). 

Proposition 3.3. {St,rit), t E [0,r], is a (time-inhomogeneous) Markov 
process with respect to P and F. Its distribution is uniquely determined by 
the SDE system (3.1)-(3.3). 
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Proof. For any h G Cb{D x {l,...,m},M) and T' € [0,r] there is a 
unique bounded classical solution v to the PDE (3.7) with terminal condition 
(3.8). By Lemma 3.2, 

E[h{ST',VT'm = E[v{T',ST',7]T')\Tt] = v{t, St,r]t) for < t < T'. 

This establishes the Markov property of {S,rj). In particular, t = gives 
E[h{ST' , rjT')] = w(5o, r/o) and so the one-dimensional marginal distributions 
of the process are unique. To show uniqueness of the finite-dimensional dis- 
tributions by induction, let /i^, . . . , /i""*"^ be arbitrary continuous bounded 
functions. For any times ti < ■ ■ ■ < tn+i < T, conditioning on Tt„ gives 



(3.10) 



E 



'n+l 

.1=1 



U^\Su,Vu)]viSt„,VtJ 



where v denotes the solution to the PDE (3.7) and (3.8) with h := /i""*"^ 
and T' ■.= tn+i- Since the right-hand side of (3.10) is determined by the 
n-dimensional distributions, the claim follows. □ 



3.3. Construction by a change of measure. At first sight, the mutual 
dependences in (3.1)~(3.3) might seem to make such models difficult to con- 
struct since we face a nonstandard SDE system where the solution (5, rj) 
also affects one part of the driving process. But the problem can be re- 
duced to the special case where TV = (7V'=J) is a standard multivariate point 
process; then rj becomes an autonomous process and S is well defined by 
(3.1). From here, the desired (t, 5')-dependent intensities for can then be 
constructed by a suitable change of measure. 

More precisely, we start with a filtered probability space (Q,^',F',P') 
carrying an r-dimensional (P',F')-Brownian motion W = (VF*)j=i , and 
a multivariate F'-adapted point process A^ = {N^-')k,j=i,...,m with constant 
(P', F')-intensity 1 for any k and j. In other words, 

/„ N^^ , k,j = 1, . . . ,m, are independent standard 

^ ' ' Poisson processes under P'. 

We assume that J-q is trivial, JF^ = JF' and F' satisfies the usual conditions. 
Then (3.2) defines an autonomous process rj. Given this process, there is a 
solution S to (3.1) under suitable assumptions on the coefficients; simple 
examples are given in the following 

Example 3.4. Provided (3.11) holds, there exists a solution {S,r]) to 
(3.1), (3.2), and the solutions X*'^'*"' to (3.5) satisfy condition (3.6) in the 
following cases [cf. Becherer (2004)]: 

(i) D = (0,00)^^, and 7,0" in (3.4) are continuous functions, depending 
only on {t,k) but not on x. In this case X*'^''^ and S can even be written 
explicitly as stochastic exponentials. 
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(ii) D = 'K'^, and F, S are Lipschitz-continuous in x, uniformly in t, for 
any k. 

For the general case, we define a probability measure P <^P' by 



where the density is chosen to give N the (P, F')-intensities (3.3); see Chap- 
ter VI.2 in Bremaud (1981). By Girsanov's theorem, VF is a local (P, F')- 
martingale whose covariance process (W) is the same under P' and P since 
it can be computed pathwise, and therefore W is also a (P, F')-Brownian mo- 
tion. Finally, if {Q.,T,¥,P) is the standard P-completion of (r2,JF',F',P), 
one can check that F satisfies the usual conditions under P. With respect to 
(P, F) we then have that VF is a Brownian motion, is a multivariate point 
process with the desired intensities (3.3) and (S, ry) solves (3.1) and (3.2). 

This change-of-measure construction extends an argument from Kusuoka 
(1999) on how to construct two point processes with mutually dependent 
intensities. Alternatively, one could infer existence (but not uniqueness) of 
a solution to the SDE system (3.1)-(3.3) from results by Jacod and Protter 
(1982); they constructed a solution for an SDE where the characteristics 
of the driving process depend on the solution process by transforming the 
problem to an SDE of ordinary type on a suitably enlarged probability space. 

4. Hedging and valuation of integrated risks. This section presents an 
approach for valuing and hedging a general class of contingent claims with 
recursive payoff structure in the model introduced in the last section, and 
shows how option values and hedging strategies for this framework can be 
obtained from solutions of suitable PDE systems of the type studied in Sec- 
tion 2. One important feature of the claims we consider is that they can be 
specified implicitly, in the sense that their own value may influence the pay- 
offs they deliver. A typical example is the pricing of a defaultable zero coupon 
bond with fractional recovery when the recovered amount depends on the 
pre-default value of the bond. Such claims lead to a fixed-point problem 
since their value depends on their payoff structure, which in turn depends 
on their value. We show how this fits naturally into the framework developed 
so far and leads, in comparison to the general setting of Duffie, Schroder and 
Skiadas (1996), to more explicit results in terms of PDEs in our setting. In 
addition to a pure pricing approach, we also offer a combination of valuing 
and hedging ideas. 
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4.1. Formulation of the problem. Our model of the financial market is 
given by the process {S,r]) from (3.1)-(3.3) with the assumptions from Sec- 
tion 3.1. We think of S as the prices of d risky assets (discounted, i.e., 
expressed in some tradable numeraire) and of rj as representing some non- 
tradable risk factors in the market. Recall that F, S and A'^-^ are in {t,x) 
and all A'^-^ are bounded. In addition, we suppose that the market price of 
risk function 

(4.1) $ := S*'(i;S*')"^r exists and is bounded on [0, T] x D x {1, . . . , m}. 

This implies that Z = £{—J ^{t, Stjijt-) dWt) is in the Hardy martingale 
space liP{P) for any p G [l,oo) and so dP = ZxdP defines an equivalent 
local martingale measure (ELMM) for S, the so-called minimal ELMM. By 
Girsanov's theorem, the dynamics of S under P are given by 

SoeD, dSt = J:{t,St,Vt-)dWt 

for a P-Brownian motion W = W + J ^{t,St, rjt-) dt, while the dynamics of 
rj and N under P are the same as under P and given by (3.2) and (3.3). 
Hence {S,rj) is also Markov under P, and 

(4.2) M''^ ■.= N''^ - J X''^{t,St)dt, k,j = l,...,m, 

are martingales in Ti^ for any p E [1, oo), under both P and P; in fact, since 
[M'^J] = iV'^J' and the intensities A''^' are bounded, a simple time change 
argument via Theorem 11.16 in Bremaud (1981) shows that [M^^]t has even 
all exponential moments. 

Remark 4.1. Except in trivial cases, there are typically many other 
ELMMs beside P so that our financial market is incomplete. Note also that 
P and P coincide for F = 0, that is, when the original probability measure 
is already a local martingale measure for S. This automatically happens if 
one starts modeling under a pricing measure, as practitioners often do. 

The financial contracts we consider are specified by functions h:Dx 
{1, . . . , m} ^ M, 5 : [0, T] X L> X {1, . . . , m} X M ^ M and f''^ : [0, T] x L> x M ^ 
M, k,j = l,...,m, with the following interpretations. [Note a slight change 
of notation in comparison to (2.22). Technically, we shall find that the func- 
tion V is given by a PDE system as in Section 2, where the interaction g 
has the form (2.22) with 6^ and f^^ depending on the argument v only via 
the kth coordinate f'^; this is sufficient for later applications and allows to 
simplify notation in the sequel. Formally, the relation to (2.22) is given by 
5^{t,x,v) = 6{t,x,k,v^) and f^^{t,x,v) = f^^ {t,x,v^) when d'', f^^ denote 
the functions from Section 2.] h describes a final payoff at time T of amount 
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h{x, k) St = X and rjx = k; 5 specifies a rate for payments made continu- 
ously in time (e.g., dividends) and f''^ describes a lump sum payment tfiat 
falls due whenever the state of r/ changes from A; to j. A typical example is 
given by a life insurance contract where rj could describe the state of health 
of the insured person. The total payoff up to time T from a triple {h,5,f^^) 
is 



(4.3) 



H = h{ST,r]T) + [ S{t,St,rjt-,v{t,St,rjt-))dt 
Jo 

+ r E f'{t,St,v{t,St,rjt.))I{,y{vt-)dN^^. 



k,j=i 



Technically, we always suppose that /i, 5 and f^^ are continuous; h is 
bounded; there is some constant K <oo such that \5{t, x, k, v)\ and \f''^t, x, v)\ 
are bounded by K{1 + \v\) for all i, x, k, j; and each 6{-, •, k, •) and all /^-^ 
are in {t,x,v) and moreover locally Lipschitz-continuous in v, uniformly 
in (t,x). By (3.2), one can rewrite the integral with respect to N in (4.3) 

as a sum EtG(o,T] Efc^j f'^K' ' ■)I{k}iilt-)I{j}{i]t)', this shows that the payoff 
involves only the random processes S and rj. Similar remarks apply at sev- 
eral points in the sequel, compare (4.4), (4.6) and (4.7), (4.13) and (4.14), 
where we use the integral notation for technical convenience. 

A closer look at (4.3) shows that H has an unusual feature. It is not (yet) 
well defined as a random variable because the extra argument v{t, St,r]t-) 
in 6 and f^^ has not yet been specified. The idea is that many intermediate 
payments depend not only on the current state of the underlying S and 
the factor rj, but also on the value of the contract itself. In life insurance, 
for instance, the payoff at death might be set equal to the current reserves; 
another example occurs with fractional recovery of a defaultable bond. How- 
ever, the function v that should give the current value of the contract still 
needs to be determined by some argument. 

A first possible idea for finding u is a pure pricing approach. We fix some 
probability measure and axiomatically define the value process of a (future) 
payment stream as the conditional expectations of all future payments un- 
der that measure. This is a common procedure in, for example, credit risk 
valuation problems. Since the modeling takes place under the valuation mea- 
sure, we may take P as that measure and thus assume that P is an ELMM 
and r = 0. By the Markov nature of the model, the valuation should be 
determined via a valuation function v, and in view of the recursive payoff 
structure we should like to write the value at time t of future payoffs as 



v{t,St,r]t) = E 



h{ST,riT) + 5{u, Su,r]u-,v{u, Su^Vu-)) du 
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,7^ m 



(4.4) 



+ / E f''{u,Su,v{u,Su,Vu-))I{k}iVu^)dNt^ Tt , 



k,j = l 



tG[o,r]. 



But since the right-hand side itself contains the function v, existence and 
uniqueness of a solution to (4.4) is a fixed-point problem; this can be viewed 
as a variant of the recursive valuation approach by Duffie, Schroder and 
Skiadas (1996). 

An alternative to a pure pricing approach is a combination of valuing 
and hedging ideas. Loosely speaking, the general goal there is to find at the 
same time a stochastic process 6 and a function v such that the trading 
strategy 9 provides a "good" hedge against the payoff H in (4.3) and the 
value process of the strategy is related to the valuation function v. To make 
this more precise, we use here as criterion for the quality of a hedge the 
concept of local risk-minimization. This goes back to Schweizer (1991); the 
idea is to find a not necessarily self-financing strategy whose final value 
coincides with the sum of all payments to be hedged, and whose cumulative 
costs over time have in a suitable sense minimal quadratic fluctuations on 
average under P. In mathematical terms, this problem can be reformulated 
as finding a decomposition of the form 



with a constant Hq and a P-martingale L that is P-orthogonal to S", where 
P is the minimal ELMM introduced above. The trading strategy which 
holds (at any time t) 9\ units of risky asset i = l,...,(i and the amount 
6^ = Hq + Lt + j^OudSu — OfSt in the numeraire used for discounting is 
then locally risk-minimizing for H under P, and one can deduce from (4.5) 
that the resulting valuation process V = Hq + / 9 dS + L is the conditional 
expectation of the payoff H under P. Up to integrability issues, this follows 
from the results in Follmer and Schweizer (1991) since S is continuous; see 
Theorem 3.5 in Schweizer (2001). When P = P, the above strategy is even 
(globally) risk-minimizing in the sense of Follmer and Sondermann (1986). 
But as with (4.4), the decomposition (4.5) cannot be obtained in a standard 
way because H involves the value function v which in turn depends via V 
on the decomposition. 

Remark 4.2. Risk-minimization is one among several hedging approaches 
for incomplete markets and (despite some drawbacks) typically leads to com- 
parably constructive solutions. Since one motivation for this section comes 
from credit risk problems where hardly any (constructive) results for hedg- 
ing complex payoffs under incompleteness seem available so far, studying 



(4.5) 
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risk-minimization is a natural first step to make. For results from another 
(nonlinear) hedging approach in a similar model but for simpler (nonrecur- 
sive) payoffs, see Becherer (2004). 

Although the reasoning behind the two approaches is quite different and 
the pure pricing approach postulates the special case P = P, both valuations 
have the same mathematical structure. Both are expectations of future pay- 
ments under the measure P, and we shall see that the valuation function v 
for both is indeed determined by the recursion formula 



v{t,St,'nt) = E 

(4.6) 



h{ST,riT) + 5{u,Su,riu-,viu,Su,'nu-))du 
+ [ f''{u,Su,viu,Su,rju-))I{k}{riu-)dNt^ 



k,j=l 



te[o,T], 



where E denotes expectation under the minimal ELMM P. To prepare for 
this, we first establish uniqueness of the corresponding valuation process 
in Lemma 4.3, whose proof is relegated to the Appendix. Existence will 
be a by-product of the main results in the next section which describe the 
corresponding valuation function v and the decomposition (4.5). 

Lemma 4.3. For any payoff triple {h,6,f^^) and probability measure 
Q € {P,P}, there exists at most one bounded semimartingale V such that 



Vt = Ec 



(4.7) 



T 



h{ST,r]T)+ 5{u,Su,r]u~,yu~)du 

H r E f'iu,Su,Vu-)I{k}{Vu-)dN^^ 
^ 1 



k,j=l 



tG[o,r]. 



4.2. Solution via interacting PDE systems. For a given payoff triple (/i, 5, f^-'), 
consider the following system of PDEs for a function v{t,x,k) on [0,T] x 
D X {1, . . . ,m}: For each fc, 

d 

= vt{t,x,k) + ^ E a^'' {t,x,k)v^i^j{t,x,k) + 5{t,x,k,v{t,x,k)) 

m 

(4.8) +Y.X''^it,x){v{t,x,j)-vit,x,k)+f''^{t,x,v{t,x,k))), 
i=i 
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€ [0,r) X D, 



and 



(4.9) 



v{T, X, k) = h{x, k) 



with {a^^)ij=i^...^d = a := SS*''. Note that (4.8) is the special case of (2.8) 
where c = 0, F = and g is given by (2.22). We want to conclude from 
Theorem 2.4 that (4.8) and (4.9) has a unique solution in Cl'^{[Q,T) x Dx 
{l,...,m},M) for any claim H of the form (4.3), and there is one point 
where we must take care. Assumption (3.6) from Section 3 guarantees non- 
explosion under P for the solution of (3.5) or (2.1) with drift F, whereas the 
present application of Theorem 2.4 requires this for the solution with drift 
0. However, this is also true and can be verified via Girsanov's theorem by 
a change to an equivalent measure under which there is no drift. 

Theorem 4.4. Under the assumptions from Section 4.1, letv € C^'^([0,T) x 
D X {1, . . . , m}, M) be the solution of (4.8), (4.9) corresponding to the triple 
{h,6,f''^). Then the payoff H in (4.3) admits a decomposition (4.5) with 



Lt= V {v{u,Su,j) -v{u,Su,k) 

+ f^^ (n, Su ,v{u,Su:r]u-W{k} iVu- ) dM^^ , 

for t G [0,T]. The process L is a martingale in TiP under both P and P, 
and {J 6dS)T is in Lp/^(P) for every p€ [l,oo). Moreover, L is a BMO- 
martingale under both P and P, and J6dS with 6 from (4.10) is in BMO 
under P. 

Proof, (i) Applying Ito's formula to v{t,St,r]t) yields 
dv{t,St,r]t) = grad^vit, St,r]t-) dSt 



Ho 

(4.10) 9t 



v{0,So,vo), 

{Ol)i=i,...,d = S'c^'^xV{t,St,r]t-) fort£[0,T) and 





m 



+ {v{t,St,3)-v{t,Suk))I{k}{m-)dNl 



t G [o,r). 
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Substituting dN^^ = dM^^ + X^^{t,St)dt, one can then use the PDE (4.8) 
to obtain 

dv{t,St,r]t) = 9t dSt + dLt - 5{t, St,i]t^,v{t, St,rjt^)) dt 

m 

(4.11) - E fKt,St,v{t,Sum-))hk}im~)'^^t' 

k,j=i 

for t G [0,r). 

The next part of the proof shows in particular that J 6dS is well defined 
on [0,T] so that (4.11) extends from [0,r) to all of [0,T] as both sides are 
a.s. continuous at time T. Admitting this for the moment, we can integrate 

(4.11) from to T and use (4.9) to conclude by comparison with (4.3) that 
the ingredients in (4.10) indeed yield the decomposition (4.5). 

(ii) To prove the desired integrability properties, we substitute dN^-' = 
dM^-' + \^^{t, St) dt in (4.11) and rearrange terms to obtain 

m 

(4. 12) Ot dSt = dvit,St,Vt) + Jtdt+ ^ K^^ dM^^ 

k,j=i 

with bounded integrands 

Jt := 5{t,St,r]t-,v{t,St,r]t-)) 

m 

+ X'^^{t,St)fHt,St,v{t,St,Vt-))I{k}{Vt-), 

k,j=l 

■= -{v{t,St,j) -v{t,St,k))I{k}{r]t_)I{k^jy 

For every p G [l,oo), M^^ is an T^^-martingale under P and P. The same 
property holds for L and for the last term in (4.12) since v, f'^^{-, ■, u(-, ■, •)) 
and K^^ are all bounded. As v and J are bounded, the right-hand side 
of (4.12) is the sum of a bounded process and an 7^^-martingale under P. 
Hence / 6 dS is well defined on [0,T] and also in 1-P{P) since it is already a 
local P-martingale. Because dP/dP has all moments under P due to (4.1), 
the asserted integrability of (/ 6dS)T under P follows via the Burkholder- 
Davis-Gundy and Holder inequalities. 

(in) Denote by U :=Y.k^j I^''^dM''^ the last term from (4.12) so that 

the local P-martingale / 6dS differs from U only by a bounded process. 
Both U and L are finite sums of stochastic integrals with respect to M^^ 
of bounded integrands, and so it is enough to prove that each (P- and P- 
martingale) M'^-^ is in BMO. But this is clear since M^^ has bounded jumps 
and (M'^'J) = / X''^{t, St) dt is bounded for both P and P. □ 
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As already mentioned, Theorem 4.4 immediately gives: 

Corollary 4.5. Under the assumptions of Theorem 4.4, the strategy 
{6^,6) holding risky assets according to 6 = {0^,...,0'^) from (4.10) and in- 
vesting 

e^:=v{0,So,i]o) + Lt+ fOudSu-etSt, te[0,T], 

JO 

into the numeraire used for discounting is locally risk-minimizing for H . If 
r = 0, this strategy is even risk-minimizing under P. 

Proof. The first assertion fohows from Theorem 3.5 in Schweizer (2001). 
The second is then clear because local risk-minimization coincides with risk- 
minimization if is a local P-martingale, that is, for T = 0. □ 

To obtain a nice representation for the valuation function u, we introduce 
the process 

Yt:=E[H\Ttl tG[0,T]. 

The structure (4.3) of the claim implies that H has all exponential moments 
under P\ this follows from the fact that h, 5, f^^ are all bounded and N^^ has 
bounded jumps and bounded intensity [see Becherer (2001), Lemma 3.4.1, 
for details]. Hence H is in LP{P) and Y is an T^^-martingale under P for 
any p G [1, oo). 

Corollary 4.6. Under the assumptions of Theorem 4.4, we have the 
representation 

v{t,St,r]t) = Yt - [ 6{u,Su,r]u-,v{u,Su,riu-))du 
Jo 

(4.13) - / J2 f'''iu,Su,viu,Su,Vu-))I{k}iVu^)dK' 

•^0 k,j=l 



E 



T 

h{ST,r]T)+ I d{u,Su,r]u-,v{u,Su,r]u-))du 



- k,j=i 

te[o,T]. 

In particular, the process Vt := v{t, St,rjt) is the unique solution to the re- 
cursive equation (4.7) for Q = P, and the function v solves (4.6). 
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Proof. By Theorem 4.4, f (0, 5o, ?7o) + / 6 dS + L is a P-martingale with 
final value H by (4.5), and so this martingale coincides with Y. Hence the 
first equality in (4.13) follows from (4.11) by integrating, and the second 
follows from the definitions of Y and H. □ 



The representation of v in (4.13) is given by the conditional expectation 
under P of future payments. This allows to view P as the pricing measure 
associated to our approach and thereby confirms (4.6). In the special case 
when the payoff functions 6 and f'^^ do not depend on the f -argument, (4.13) 
can also be used to compute v by Monte Carlo methods and thus provides 
an alternative to the numerical solution of the PDE system (4.8) and (4.9). 

Remark 4.7. To motivate how the PDE (4.8) arises, let us recall from 
Follmer and Schweizer (1991) that (4.5) is basically the Galtchouk-Kunita- 
Watanabe decomposition of H under P. For the martingale K = S[i?|F], 
the structure of H in (4.3) yields 

Yt= [ 6{u,Su,'nu^,v{u,Su,'nu~))du 



(4.14) 



+ / Y.f'''(^^Su,viu,Su,r]u-))I{k}{Vu-)dN^^ 



k^3 



+ E 



T 



h{ST,'nT) + 6{u,Su,i]u-,v{u,Su,r]u^))du 



+ 



T 



E 



f''^{u, Su, v{u, Su, Vu~))Iik}{Vu~) dN^^ 



Since {S, rj) is a Markov process under P, the last term is a function of 
{t,St,r]t) only, and we call it v{t,St,r]t). Assuming that v is sufficiently 
smooth, we can apply Ito's formula, and since y is a P-martingale, all drift 
terms on the right-hand side must add to 0. Writing out the zero drift condi- 
tion then produces the PDE system (4.8). However, this is only a heuristic 
argument because we have given no reason why v should be sufficiently 
smooth, nor why the last term in (4.14) should be equal to v{t, St,rit). The 
latter is a genuine problem since v appears at many places in (4.14), and 
this explains why a proof requires a fixed-point argument. 



5. Examples, applications, extensions. This section outlines some possi- 
ble applications of our general model. We briefly sketch some links to insur- 
ance risk problems and dwell in more detail on issues related to credit risk. 
There we first explain how the contributions made here fit into an overall 
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perspective, and then show how one can extend the results to incorporate 
further important aspects in the context of credit risk: hedging by using 
default-related assets, and stock prices that crash at default. 

5.1. A hybrid reduced form model for credit risk. In the context of credit 
risk modeling, it is natural to have a finite-state process ij representing the 
rating (or default) state of i entities under consideration; we think of rjt G 
{AAA, . . . , C, D (default)}^. As usual in reduced form models, the time of 
default is defined by a jump of a point process which drives rj. The intensities 
for rating changes may depend on the stochastic evolution of the Ito process 
S. If the SDE coefficients for S do not depend on r/, then S" is a diffusion 
process, is a Cox process and 77 is a conditional Markov chain. Then the 
model falls into the Cox process framework of Lando (1998) which in fact 
inspired our model. Alternatively, it can be viewed as a (Markovian) reduced 
form model with state variable process S; see Bielecki and Rutkowski (2002) 
for more explanations and many more references. 

Remark 5.1. One can ask whether tradable asset prices should be af- 
fected by the rating from agencies, or rather by some nonobservable solvency 
state of the firms. We do not address this question here, but use the rat- 
ing process only as an example for an observable finite-state process which 
reflects credit risk. One reason to consider the official rating may be the ex- 
istence of financial products whose payoffs are linked to a rating by agencies. 

In comparison to existing credit risk literature, this article makes several 
contributions. Because of the mutual dependence between assets S and rat- 
ing process r], the processes and r], conditional on the evolution of S, 
are in general not a Poisson process respectively a Markov chain with in- 
tensity matrix {X^^ {t, St))kj- In this regard, the modeling goes beyond the 
Cox process setting of Lando (1998). Moreover, the model is hybrid in the 
sense that it incorporates prices S of liquidly traded assets to which the 
default risk is related, with default intensities being a function of (t,St). 
Such ideas go back to Madan and Unal (1998). If the default intensities of 
certain firms depend, for instance, on the stock indices of related industries, 
these indices could be used for hedging and must be taken into account 
for pricing. Furthermore, we let the payoffs of contingent claims depend on 
tradable assets and on the default process rj, and we can to some extent also 
deal with products like convertibles which are related to the stock price of 
a firm itself; see Section 5.1.2. Yet another contribution is that we consider 
not only a pure pricing approach under an a priori pricing measure as is usu- 
ally done in the credit risk literature. We also offer an alternative combined 
hedging and valuation approach with the aim of minimizing hedging costs. 
Note that the model should then be specified under the probability measure 
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under which hedging costs are to be minimized. Finahy, the solutions to the 
pricing and hedging problems are described in terms of PDEs under precise 
conditions; we do not just assume "sufficient regularity" as is often done in 
the literature. 

To exemplify some features of our model, we first consider a single de- 
faultable entity with 

r]t G {no default (n), default (O)} 

and default being an absorbing state. If we choose as discounting numeraire 
the zero coupon bond with maturity T, prices are T-forward prices and the 
payoff of a defaultable zero coupon bond is /i(-S't, ??t) := I{n}{'nT)- Fractional 
recovery in the case of default by some fraction ii € [0, 1) from treasury or 
from pre-default (market) value is modeled by choosing f"^{t,x,v) := R or 
:= Rv in the payoff specifications (4.3). 
For basket products or counter-party credit risk, one has to deal with 
multiple default risks and modeling the dependences becomes a crucial task; 
see Chapter 10 in Schonbucher (2003). In our setting, a natural choice for ij 
is the joint rating state of the I firms under consideration, that is, 

{r]l)i=i,...,i = {AAA, AA, . . . , C, D}^ I e N. 

The intensity for a change from the current joint rating state r^t- to another 
state j is then 

X'^^it, x)\k=^,_^,=s, for j G {AAA, . . . , C, D}^ 

and can be specified as a function of time t, current asset prices St and the 
current overall rating rjt- itself. This modeling gives a flexible Markovian 
framework which includes several more specific parametrizations suggested 
so far. It permits simultaneous defaults of different firms, default correlation 
because of joint dependence of individual defaults on the common (factor) 
process 5", as well as sudden jumps in the default intensities of individual 
firms when other related firms default. For an implementation, one has to 
decide on a reduced version with as many or few parameters as can be 
reasonably fitted from available data. 

For a more concrete example, consider I default risks and for r/ the sim- 
plified state space {n, D}^ with no further ratings. We exclude simultane- 
ous defaults and make default an absorbing state for each risk by setting 
\^^{t,x) := when J2i=i^{k''j^j'} > 1 or when {k^,f) = (f,n) for some i. 
Then we are left to choose X^^{t,x) for those k ^ j with {k^,f) = {n,d) for 
exactly one i = i* and /c* = f elsewhere. As pointed out in Davis and Lo 
(2001) and Jarrow and Yu (2001), defaults of some firms may trigger jumps 
in the default intensities of other firms, and a simple way to capture this 
phenomenon is to take 

A''^(t,j;) = AaSti^{n('=')) with A > and a > 1 



REACTION-DIFFUSION SYSTEMS AND HEDGING 



27 



for those k^j with = {n,d) for exactly one i = i* and A:* = f for 

all other i. The basic default intensity of any firm is then A and increases 
by a factor a > 1 whenever another firm has defaulted. Modeling default 
dependences between £ firms is thus reduced to two parameters. To obtain 
higher default intensities when stock indices are low, practitioners have sug- 
gested to model the default intensity as a function of related tradable assets; 
see Arvanitis and Gregory (2001) or Davis and Lischka (2002). Pushing the 
example a bit further, one could incorporate such effects in our model by 
taking A = A(t, St) as a function of a single (d = 1) tradable index S. A next 
step could be a factor model with different groups of firms, and so on. 

Instead of considering now one particular application and writing down 
the resulting PDE system for that problem, we explain in the following two 
subsections how to extend our framework to address additional aspects of 
pricing that seem relevant in the context of credit risk. 

5.1.1. Perfect hedging in a completed market. Up to now, we have as- 
sumed that S represents the price processes of all tradable assets which are 
available for dynamic hedging. But for some applications, it is natural to 
suppose that there are additional defaultable securities which are liquidly 
traded, and one feels it should be possible to replicate a given defaultable 
contingent claim if there are enough tradable securities which are sensitive 
to the default event. We now show how to make this intuition precise. 

To explain the idea in the simplest setting, we consider a single firm which 
can only be in default or not in default. So we take the model from Section 3.1 
for the process {S, rj) with St £ D = (0, oo) and r]t G {no default (n), default (O)}. 
The default state is absorbing so that A®"(t,x) := 0. If r/o = n so that the 
firm is not in default at time t = 0, the default time is 

r := inf{t € [0,r]|?7t = c)} = inf{t G [0,T]|AiV4"V 0}. 

We suppose that the drift of S vanishes, that is, ^{t, x, k) := for all t, x, k, 
and we think of P as an a priori pricing measure. But in contrast to Sec- 
tion 4.1, we assume here that there is in addition to S a further tradable 
security which is sensitive to the default event. For concreteness, we choose 
as discounting numeraire the (nondefaultable) zero coupon bond with matu- 
rity T and suppose that the defaultable zero coupon bond with maturity T 
and recovery from treasury is tradable. By Theorem 4.4 and Corollary 4.6, 
the dynamics of its market price process Yt := E[I^.,-^x} + RI{T<T}\^t] with 
i? € [0, 1) are given by 

dYt = dv{t, St,vt) = v.,{t, St,vt-) dSt + {v{t, St, 0) - v{t, St, n)) dM^' 

and are described by the solution v^{t,x) = v{t,x, k) to the PDE system 

^ = Vt{t,x) + \x'^a'^{t,x,k)vl^{t,x) + I{k=n}>^''^{t,x){v^{t,x) -v'^{t,x)) 
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with {t,x) G [0,r) X (0,00), k = n,d and boundary conditions v^(T,x) = 
-^{n}(^) + RI{d}ik). In the present setting, v^{t,x) = R, and the payoff of a 
defaultable claim of the form (4.3) is 

H = h{ST, nr)+j^ S{t, St,vt~,vit, St,vt-))dt + I^r<T}r\r, Sr, v{t, 5., n)). 

Again by Theorem 4.4 and Corollary 4.6, the dynamics of Yt := E[H\J^t], 
te [0,r], are 

dYt = v^{t, St,vt-) dSt + {v{t, St,d) - v{t, St, n) + r\t, St,v{t, St, n))) dMf, 
where v^{t,x) =v{t,x,k) denotes the solution to the PDE system 
= Vt{t,x) + ^x'^a'^{t,x,k)v^^{t,x) 

+ 6it,x,k) + (A;) {t, x) {v^ (t, x) - r;" {t, x) + (t, x, x))) 

with {t,x) G [0,r) X (0,00), k = n,d and boundary conditions v^(T,x) = 
h{x,k). From the dynamics of Y and Y, we conclude that 



Vx{t,St,r]t-) dS, 



dYt = [vx{t,St,rit-) 

v{t, St,D) - v{t, St, n) + r\t, St,v{t, St, n)) 
v{t,St,d)-v{t,St,n) 

f v{t,St,d)-v{t,St,n) + r\t,St,v{t,St,n)) \ - 
\ v{t,St,d)-v{t,St,n) J 

Hence, the claim H = Yt can be replicated by self-financing dynamic trading 
in the risky securities S and Y , and the hedging strategy is described in terms 
of the solutions v and v to some PDEs. By absence of arbitrage, the value 
process of H is then given by the wealth process Y of the replicating strat- 
egy. This can be viewed as a concrete example (describing the hedging strat- 
egy by a PDE) for the more abstract and general representation theorems 
by Belanger, Shreve and Wong (2004) and Blanchet-Scalliet and Jeanblanc 
(2004). For other results on the hedging of replicable defaultable claims, see 
the recent article by Jeanblanc and Rutkowski (2003). 

5.1.2. Stock prices with downward jumps at default. So far, we have mod- 
eled by S the price processes of tradable assets and since S is continuous, 
we interpret them as stock indices of industries related to the considered 
defaultable entities. If we consider the tradable stock of a firm itself, things 
change because we expect this price to jump downwards when the firm de- 
faults on some of its payment obligations. Since products like convertible 
bonds involve both the stock price and the default state of a firm, this case 
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is very relevant for applications. This section shows how to extend our frame- 
work to deal with such situations and obtain similar PDE results as before. 
We explain the ideas in the simplest setting with a single firm which can be 
either in default or not, and we restrict ourselves to a pure pricing approach 
without considering hedging issues here. 

We start with the model of Section 3.1 for (5, r/) with r]t G {n,l)} and 
St € D = (0, oo). The basic idea is to view S as the pre-default value of the 
stock. The firm is not in default at time t = and default is an absorbing 
state so that r/o := n and X^"{t, x) := 0. Default happens at the stopping time 

(5.1) T :=mi{t e [0,r]|77t =0} =inf{t G [0,T]\AN^' + 0}, 

and the events of default or no default up to time t are given by {t > r} = 
= f } and {t < r} = {rn = n}. For a fractional recovery constant TZ G [0, 1], 
the firm^s stock price S is 

(5.2) St := Sti{t<r} + nsti{t>r} = St-{i- n)Sti{t>r}, t G [o,r]. 

Hence the stock price drops to a fraction TZ of its pre-default value when 
default happens. By (5.2) we have dS = S^{dS - (1 - TZ)I{„y{r]) dN"^) and 
by the SDE (3.4) for S we conclude 

dSt = St-{{j{t,St,Vt-) - (1 - n)I{,}{rit.)X''^{t,St))dt 

+ a{t, St, m-) dWt - (1 - 7e)/|„}(7?t_) dM,"°) 

with the martingale M"' as in (4.2); this uses that St = St- = St- on 
{•qt- = n}. Since we want P to be a pricing measure, S should be a lo- 
cal P-martingale, and if we start from a volatility function a and a default 
intensity function A"^, this requires to take as drift 

7(^,x,A:):=(l-7^)/{„}(A:)A"'(^,x). 

With this choice, the SDE for S can be written without reference to S as 

dSt = St-{a{t, St-,rjt-) dWt - (1 - 7^)/|„|(r?t_) dMf ), 

with a{t,x,n) := cr{t,x,n) and a{t,x,D) := a{t,x/TZ,d) if 7?, > or a{t,x,d) 
arbitrary if 7^ = 0. It is clear that S is then a local (P, F)-martingale. To 
derive the pricing PDE with P as pricing measure, we start for TZ £ [0, 1] by 
considering claims of the form 

H = h{ST)I{r>T} + Ht, St, v{t, Sr-)) dt + /(r, Sr-,V{T, Sr-))I{r<T} 

with sufficiently regular payoff functions h, 6 and / (precise conditions are 
given below). As before, v denotes a valuation function which gives the 
value of outstanding payments from the claim, given that default has not 
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yet happened, and which stih has to be determined. Using definition (5.1) 
of r and the equality St- = St^ = St on JO, r] , one can rewrite the claim as 

H = h{ST,VT) + [ 6{t,St,r]t_,v{t,St,rjt-))dt 
Jo 

+ r I{„^{r]t-)r'{t,Suv{t,Sun))dN:t' 
Jo 

with h{x,k) =h{x)I^i,^ny,6{t,x,k,v) =6{t,x,v)I{k=n}^ r^{t,x,v) = f{t,x,v) 
and v{t,x,k) = v{t,x)I^f^=n}- This represents in a form like (4.3) which 
refers to S and r] instead of S and r. Now suppose that h, 5 and /"'' sat- 
isfy the regularity conditions stated after (4.3) and recall from Lemma 4.3 
(with Q = P) that there is at most one valuation process for a claim with a 
recursive payoff structure like H. One can then show that the price E[H] at 
time equals v{Q,So,n), where v is (see Theorem 2.4) the unique bounded 
solution to the PDE system 

^ = vt{t,x,n) + (1 -7^)A"'(^,a;)3;^;a;(^,x,n) 

+ ha^it^x^n^x^VxxitiX^n) + 5(t, x, n, u(t, x, n)) 

(5.3) 

+ A"°(t, x){v{t, X, 0) - v{t, X, n) + r°(t, X, v{t, x, n))), 

= vt{t,x,d) + ^a'^{t,x,d)x'^Vxx{t,x,li), 

for {t,x) G [0,T) X (0,oo), with boundary conditions v{T,x,k) = h{x,k) for 
k = n,d. Since v{t,x,d) = is trivial, we are left with one (simplified) PDE 
for t;(t,x,n). For fractional recovery of the claim in default, that is, /"^(t,x,f) = 
Rv with recovery constant R S [0, 1] (for the claim), the last term in (5.3) 
simplifies further to — A"''(i, a;)(l — R)v{t,x,n), exhibiting the well-known 
structure "default intensity times loss fraction." The proof of the above val- 
uation result goes as for Theorem 4.4; one applies Ito's formula to v{t, St-,r]t) 
to show that 

Yt = v{t,St,r]t) + I 5{u,Su,'nu-,v{t,Su,riu-))du 
Jo 

+ f /{„}(77„_)/"'(^x, v{t, Su, n)) dTVr 
Jo 

is a martingale with final value Yx = H so that E[H] =10 = ^^(0, 5*0, %)• 

In the case when TZ is strictly positive, the stock price at default cannot 
drop to zero. By (5.2), we then have not only St = St~ on [0,r]l, but also 
St = St- /TZ on Jt, TJ. This permits to consider more general claims of the 
form 

H = h{ST,r]T) + [ S{t,St-,rjt-,v{t,St-,'nt-))dt 
Jo 
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+ I{r<T} f St- , v{t, Sr-,Vr~)) 

= h{ST,r]T) + [ S{t,St,'nt-,v{t,St,r]t-))dt 
Jo 

+ r I^n}{r]t~)r\t,St,v{t,St,Vt-))dN^\ 
Jo 

where the functions with argument 5 or 5 are related by h{x, k) = Is^„y{k)h{x, n) 
/{J,} {k)h{x/n, D), 5{t, X, k, v) = I{„} {k)6{t, X, n, v) + I[^y{k)5it, x/7^, c), v) , /"^ (t, x 
f{t, X, v) and v{t, x, k) = v(t, x, k)I^^^{k) + v{t, x/TZ,d)I^jiy{k). If h, 6 and 
satisfy the assumptions after (4.3), arguments analogous to those above yield 
the price E[H] = v{0, So,n) for v now solving 

= vt{t, X, n) + (1 - 7^)A''°(^, x)xv^{t, x, n) 

+ ^a'^{t,x,n)x'^Vxx{t,x,n) + 6{t,x,n,v{t,x,n)) 

+ A"°(t, x){v{t, X, X)) - v{t, X, n) + x, v{t, x, n))), 

= Vt{t, X, O) + |cr^(t, x,d)x'^Vxx{t, x,d) + 6{t,x,d, v{t, x,d)), 

for {t,x) e [0,T) X (0,oo), with boundary conditions v{T,x,k) = h{x,k) for 
k = n,d. This system is in general more complicated because v{t,x,d) will 
no longer vanish if the claim H comprises payments after default time r. 

5.2. Insurance risk. A second major area for examples and applications 
is the treatment of risks at the interface of finance and insurance. This 
has attracted much attention recently, and we briefly mention a few links to 
the present work. We concentrate on valuing insurance-related products and 
leave aside other important topics such as portfolio optimization for insurers 
or ruin probabilities in the presence of financial markets. 

The basic idea is to study a model containing both financial and actu- 
arial components where S describes as before the tradable assets in the 
financial market, while r] is now related to the evolution of some insur- 
ance contract. For a concrete example, think of ijt as the state of health 
at time t of a policy holder; this may include death if we consider a life 
insurance policy. The functions h, 5 and / specify the contract's final pay- 
ment h at expiration, continuous payments (e.g., premiums) at rate 6 and 
lump sum payments / that occur on passing from one state to another. 
One example of a recursive payoff here arises if death benefits are a frac- 
tion of the reserve; see Ramlau-Hansen (1990). In general, the goal is to 
determine a value for the contract given by {h, 6, f) and possibly also to 
hedge it. 

While such products per se are not new, recent developments have started 
to put emphasis on their valuation by market-based methods that go beyond 
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traditional actuarial approaches. Examples include risk-minimization or in- 
difference pricing [see M0ller (2001, 2003)] or computation of market-based 
reserves [see Steffensen (2000)]. In comparison with existing work, we of- 
fer here two contributions. We abandon the assumption of independence 
between the financial and actuarial risk factors imposed in M0ller (2001, 
2003), and give a rigorous construction for the model with mutual depen- 
dences between S and rj. Such dependences can be important for products 
in non-life insurance; a catastrophic insurance event might, for instance, 
affect stock indices of related industries. Moreover, we provide precise con- 
ditions and existence results for PDEs and fixed points. This contributes to 
the study of related actuarial applications in Steffensen (2000) and gives a 
sound mathematical basis to previous intuitive ideas. It will be interesting 
to see how the methods developed here can be used further, for instance to 
prove verification results for other problems in the area of insurance. 

APPENDIX 

This appendix contains some proofs that were omitted from the main 
body of the article. 

Proof of Lemma 3.2. This is similar to the first step in the proof 
of Theorem 4.4. We apply Ito's formula to v{t,St,rjt), substitute dN^^ = 
dM^^ + \^^{t,St)dt and use the PDE (3.7) to conclude that the drift term 
vanishes on [0,T'). This yields 

dv{t, Sum) = ((grad, vf'T.){t, St, r,t-) dWt 



It', 



+ J2x''^{t,x)iv{t,St,j)-v{t,St,k))I{,}{7jt-)dM^ 
i=i 

tG [o,r'). 

Since v and A'^-^ are bounded and the martingales M^^ are in TiP{P) for any 
p € [1, oo) by the remark following (4.2), all stochastic integrals with respect 
to M^^ are in 'W{P), and so this must hold for the one with respect to W 
as well. In particular, the above equation extends to [0,T'] and the process 
v{t,St,rjt), t G [0,T'], is a martingale. □ 

Proof of Lemma 4.3. Suppose and V'^ are bounded semimartin- 
gales which both satisfy (4.7). By using the martingale compensator X^^{t, St) dt 
for the jump process N^^ [which is the same under P and P, see (4.2)], the 
recursive representation (4.7) of can be rewritten as 



(A.l) * 



T 



h{ST,VT)+ g{u,Su,Vu-,V^-)du 



Tt 

iG [0,r],i = l,2, 
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with g{t,x,k,y) := 5{t,x,k,y) + J2j^k^^^{t^^)f''^it^x^y)- Since all •) 
and f^^ are locally Lipschitz in y, uniformly in {t,x), g is also locally Lips- 
chitz in y, uniformly in {t,x,k). If < oo is an upper bound for \V^\, \V'^\, 
the local Lipschitz property gives an L = L{K) < oo such that \g{t, x, k,y^) — 
g{t,x,k,y'^)\ < L\y^ — y'^l for any < K, i = 1,2, and all {t,x,k). Similarly 
as for Proposition 2.1, it follows from (A.l) that we have for any (3 G (0,oo) 



-l3iT-t)>yl 



<e-/5(^-*)i?, 
L 



\g{u, Su,r]u-,V^-) - g{u, Su,r]u-,V^-)\du 



sup e-^(^~")L|l/ 

«G[0,T] 



< 



sup e-^(^-")|K' 
«e[o,T] 



So if the RCLL semimartingales and V"^ were not indistinguishable, 
considering the supremum of the above inequality over t € [0,T] would lead 
to a contradiction for jS > L. □ 



The second part of the above proof varies an argument from Duffie and Epst 
(1992) (Appendix with Skiadas). Their result cannot be applied directly be- 
cause the function g above is only locally Lipschitz. And the exponential 
weighting permits to shorten the proof. 

Acknowledgments. The first author thanks Fredi Troltzsch, TU Berlin, 
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